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Abstract 

We study the feasibility of a sterile neutrino search at the China Advanced Research Reactor by 

measuring v e survival probability with a baseline of less than 15 m. Both hydrogen and deuteron 

have been considered as potential target. The sensitivity to sterile-to-regular neutrino mixing 

is investigated under the "3 (active) +1 (sterile)" framework. We find that the mixing parameter 

sin 2 (2#i4) can be severely constrained by such measurement if the mass square difference Am 2 4 is 

of the order of ~eV 2 . 
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I. INTRODUCTION 



Neutrino flavor mixing and oscillation, a direct consequence of non-zero neutrino masses, 
are well established by experimental data |l|. Most of the experimental data to date, includ- 



ing the recent discovery of 6*i 3 



4j at nuclear reactors, can be described under the 3-flavor 



mixing framework [1], which fits naturally into the Standard Model of particle physics. 

Beyond the standard 3-flavor model, sterile neutrinos are postulated as a special type of 
(heavy) neutrinos that do not interact electromagnetically, weakly, or strongly, hence the 
name "sterile". Since 90's in the last century, there have been several neutrino oscillation 
experiments 5 ] ] jv ] js ] which seemed to detect anomalies beyond the 3-flavor mixing. Quite 
recently, after a re-evaluation of the nuclear reactor flux prediction [9], a global deficit (2- 
3 sigma level) is emerging in the measured flux from all short baseline reactor neutrino 



experiments 



101 ] . Global fits under "3 (active) +1 (sterile)" (or "3+1" in short) framework 



111- 



favor a sizable sin 29^ ~ 0.1 and a mass splitting Amj 4 ranging from 1 to 5 eV 

Although somewhat non-standard, sterile neutrinos appear to be the simplest explanation 

of existing experimental anomalies. They are also candidates for warm or cold dark matter 

under many theoretical models 12 ]. 

Numbers of experiments are underway worldwide to search for normal-to-sterile neutrino 

oscillations [131 ] . including several projects at nuclear reactors. In this paper, we evaluate 

the feasibility of carrying out a short baseline neutrino experiment using the upcoming 

state-of-art China Advanced Research Reactor (CARR) reactor 14j with a thermal power 

of ~ 60 MW th . 



II. STERILE NEUTRINO SEARCH AT CARR 



As will be demonstrated below, to search for eV-scale sterile neutrinos via reactor neu- 
trinos (a few MeV in energy) disappearance, short baseline (< 15 m) is needed to have 
sufficient sensitivity. Due to safety regulations, it is nearly impossible to place the detector 



so close to commercial reactors (~GW^ 
smearing effects to the oscillation signal 



. In addition, the core size (~meter) introduces 
151 ] . Compact research reactors (~0.1 GW t /J, on 
the other hand, are more advantageous in these regards (with a cost of lower neutrino flux 
from the core). 
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A. CARR experimental site 



China Advanced Research Reactor (CARR), constructed at the China Institute of Atomic 
Energy in Beijing China, is a tank-in-pool, inverse neutron trap type, light water cooled, 



heavy water reflected, multi-purposes research reactor [14j. The reactor body is immersed in 
a water pool with 16 m in depth and the core is located 12 m below the pool water surface. 
The reactor core is about 0.85 m in height and 0.40 m in diameter. CARR takes U3Si2-Al as 
the fuel meat, with a 20% enrichment of 235 U in weight. With a thermal power of 60 MW, 
the maximum output thermal neutron flux is 1.0 x 10 15 n/cm 2 /s. 

CARR was designed as a general thermal neutron facility for material and biological re- 
searches, as well as for isotopic production/enrichment. At present, there are 9 horizontal 
and 21 vertical beam lines coupled to user equipments. On the ground level of the exper- 
imental hall (where horizontal beam lines are), the center of the reactor core is 120 cm 
above the floor. The outer diameter of the concrete shielding structure is 5.0 m. With this 
geometry, the closest radial location for a neutrino detector is about 7 m from the core. 

B. Neutrino flux and spectrum at CARR 

Nuclear reactor is a very intense source of neutrinos. Pure electron antineutrinos z/ e s are 
produced via /3-decay of fission fragments. For a 1 GW t h reactor, there are approximately 
2 x 10 20 z/ e s emitted per second. 

Fission nuclei are dominated by 4 isotopes, 235 U, 238 U, 239 Pu, and 241 Pu, while the other 
isotopes contributes less than 0.1%. Direct theoretical calculation of the neutrino flux and 
energy spectrum bare large uncertainties (at the level of 10%) 9], primarily due to incom- 
plete information from nuclear databases. On the other hand, for 235 U, 239 Pu, 241 Pu, v e 
energy spectra have been derived from measured electron spectra at ILL, with an average 
uncertainties less than 2%, mainly originated from the uncertainty due to the conversion 
from electron to neutrino spectra 16|. For 238 U (fast neutron-induced fissions), only theo- 
retical calculations exist at present jjj, [l^ . The isotopic concentration in fuel evolves with 
reactor operation time, so does the fission rate of each isotope. To predict neutrino flux 
at a given time, commerical reactor in particular, detailed core simulation is needed, which 
carries its own uncertainty. For research reactors like CARR, however, 235 U enrichment is 
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much higher than that of commerical reactors, fissions of which dominate the total fission 
rate. For simplicity, we shall ignore burnup effects and assume a pure 235 U neutrino spec- 
trum for CARR in the remainder of this work. The difference in shape between a real and 
a pure 235 U will be considered as a shape uncertainty (bin-to-bin). 



In this work, we adopted a simple parameterization 18) for 235 U neutrino energy spec- 
trum, 

f(E ) = e a870 -°- 160S "- a09iaE v (1) 
in units of P e /(MeV • fission), also shown in Fig. [IJ 

CARR neutrino flux 

x10 ,s 




FIG. 1: Emitted neutrino flux (/s/MeV) at CARR assuming pure 235 TJ fissions with the 
parametrization given in [181 ] . 



The energy release per fission for 235 U is 201. 7± 0.6 MeV in [19| and 201.92 ±0.46 MeV in 



20j. We take an average with E re = 201.8 MeV. For CARR (P th = 60 MW), the expected 



neutrino spectrum emitted from the core per unit time is: 



F(E V ) = ^ ■ f(E v ) . (2) 

Hire. 



III. v e DETECTION 

Neutrino-target interaction cross section is low, typically of the order of 10 -44 cm 2 . Back- 
ground suppression is a key consideration of such experiments. On the other hand, short 
baseline requirement and space constraint dictate that the neutrino detector has to be placed 



4 



above ground without large shielding structure. To supress background, we investigated 
three detection techniques, all with timing coincidence characters. 



A. Inverse-/? decay with liquid scintillator 

The classical method for detecting reactor v e is the so-called inverse-/? decay (IBD), 

v e + P — > e + + n . (3) 

Liquid scintillator (LS), usually with >10% hydrogen in mass (11% assumed in this study), 
is commonly used both as the target and detector for this reaction. The positron loses 
energy (kinetic energy and the 1.022 MeV annihilation energy) immediately, emitting prompt 
scintillation lights. The neutron will undergo thermalization collisions with hydrogen, and 
eventually get captured, emitting gamma rays which are converted into delayed scintillation 
lights. The neutrino energy can be reconstructed via E e + = E 1/ — (M n — M p ) ~ E v — 1.8 MeV, 
where E e + and E v are the kinetic energy of the positron and neutrino, respectively. The 
detected prompt energy E p , on the other hand, contains both the positron kinetic energy as 
well as the annihilation energy, i.e. E p = E e + + 1.022 MeV. To enhance neutron detection 
efficiency and suppress background, most experiments adopts Gadolidium-doped LS (GdLS) 
- the ~8 MeV n-Gd capture gamma ray can be used as a clean neutron tag. 

Taking into account higher order electroweak corrections, the cross section of the IBD is 



given m 



2l|. Folding it with reactor neutrino spectrum (Eqn.[2]), detected neutrino spectrum 



(without oscillation) can be written as 

N 

N no - osc {E u ) = —J^ e (E u )E(E u ) x a IBD {E v ) x T , (4) 

where N p is the number of target protons, e(E u ) is the detector efficiency, and T is the 
duration of the measurement. The resulting neutrino spectrum is shown in Fig. |2} where the 
1.8 MeV reactor threshold of the IBD is manifest in the curve. To set the scale, the average 
detected neutrino rate (assuming 100% detection efficiency) is about 7000/MW/ton/year at 
7 m. The oscillated neutrino rate and spectrum will deviate from this spectrum for a given 
set of oscillation parameters, as illustrated in Fig. |2j 
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Detected neutrino spectrum for IBD 
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FIG. 2: Measured neutrino spectrum of IBD without and with oscillation. Uncertainties are 
statistical only. 

B. IBD with light water 



A common background for IBD detection with LS is the fast neutron background. The 
recoiling protons created by a fast neutron will scintillate and mimic the prompt energy from 
the IBDs before the neutron get captured. Since CARR is a surface facility, fast neutrons 
background induced by comic muons and from the reactor itself may pose serious challenge 
to the experiment. 

To mitigate this, if the target is water instead of scintillator, recoiling protons will not be 
able to make Cerenkov lights, therefore get rejected. The technology of Gadolinium-doping 
in water has been proposed in 2004 |22j and under active development [23|, so it is possible 
to maintain the clean neutron capture tag. However to use this approach in reactor neutrino 
experiment there are two obvious issues. 

1. The amount of Cerenkov photons is much smaller compared to the scintillating lights. 
As a result, water detector will have much worse intrinsic energy resolution compared 
to that of LS. We would have to increase the photocathode coverage to get a reasonable 
amount of photoelectrons (PEs). 

2. The Cerenkov threshold for water is about 289 keV (positron kinetic energy). The 
(nonlinear) reconstruction from the visible Cerenkov lights to true neutrino energy 
will be quite different from the LS case. 
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A realistic estimate of the Cerenkov light yield is needed to address the first concern. Light 

n 

yields for large water detectors are summarized in [l[ with a range between 3-9 PE/MeV. 
Small detectors with less light attenuation and more reflections could end up with more 
lights. To test this, we performed a bench test using tagged cosmic ray impinging on an 
acrylic ball with 5 cm diameter containing pure water. The Cerenkov lights are viewed by 
four Hamamatsu R7725 photomultipliers close-by. The measured PE spectrum is shown 
Fig. |3j The minimum ionization bump (~8 MeV) is located at ~33 PE. If we extrapolate 
this result (photocathode coverage of ~12%) to an experiment with >50% coverage, a light 
yield of 16 PE/MeV energy deposition is attainable. The energy resolution is still quite low 
compared to the LS, but it appears that a 30% resolution is achievable at 1 MeV. 



Spectrum Using Water 
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FIG. 3: Measured muon photoelectron spectrum with a prototype water detector. The solid angle 
coverage by PMT cathode is ~12%. 



The second concern can only be addressed through comprehensive calibration. Not only 
one needs a careful calibration with gamma sources, but also electron sources (beta or con- 
version electrons) to establish the relationship between the gamma energy peak to that of an 
electron at the same energy. The elaborated calibration programs developed at SuperK js4 1 



and SNO 



25] provide invaluable guidance in this regard. 



C. z/-D with heavy water 

Room gamma background, when in accidental coincidence with random (cosmic or reac- 
tor) neutron capture signals, can form IBD-like background. If such background is significant 
compared to the IBD signals, it may become a serious issue, although in principle such a 
background can be statistically subtracted. 

To further suppress background, we have considered heavy water (D 2 0) as a potential 
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target. In 70 's of the last century, Reines et al. pioneered v-D measurements at nuclear 
reactors 



261 ] . Reactor z/ e s are detected via z/ e -deuteron charge current scattering: 

u e + D^n + n + e + (5) 

with the total energy of positron given by 

E e+ ~E U - 4MeV (6) 

The detection signal now becomes a triple coincidence between the prompt positron signal 
and two delayed neutron capture signals. Neutron captures on deuteron will give a single 
6.25 MeV gamma ray. To avoid energy leakage for this high energy gamma, one could dope 
the heavy water with salt (NaCl) so that neutron capture on Clorine gives a total energy of 



8.6 MeV, distributed in 2 or 3 gamma rays 



The signals produced by fast neutron will 



be singles, the same as in a water detector. Accidental backgrounds will be highly suppressed 
by the triple-coincidence requirement. It should be noted that v e can also scatter off from 
D via neutral current channel, v e + D — > n + p + u e , but with no coincidence signature. 



A tabulated charge-current z/ e D cross section can be found in [28|. The measured spec- 
trum with no oscillations is now 

Nn 

N no - osc (E u ) = -JL. t (E u )F{E u ) x a„- D (E v ) x T . (7) 

With a lower cross section and higher energy threshold, the total number of detected neutrino 
events at 7 m is about 92/year/MW/ton (100% detection efficiency), 2 orders of magnitude 
less than that of IBD. Illustrated in Fig. H] is a comparison of a non-oscillation and oscil- 
lated spectrum, with realstic statistical uncertainty. One should emphasize that the main 
advantage of D 2 is its highly unique v-D charge current signature which may lead to a 
background-free measurement. 



IV. BACKGROUND ESTIMATION 



Direct determination of reactor-on background in the CARR experimental hall has yet 
to be performed, awaiting for the operation of the reactor. In this section, background 
estimation will be given based on current best knowledge at CARR and projections from 
other experiments, serving as a rough guidance to the design of the experiment. 
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Detected neutrino spectrum for vD 
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FIG. 4: Measured z^-D charge current spectrum without and with oscillation. Error bars are 
statistical only for a 1-ton target at 7 m operating for a year. 

A. Trigger rate 

The raw trigger rate for all three type of detectors will be dominated by room gamma 
background. We have not made dedicated gamma spectrum measurement at CARR, but 
earlier commissioning run indicated that the dose rate is less than 3 /xSv/h, which translate 
to an upper limit of gamma flux of 150 Hz/cm 2 , or <~7.5 MHz for a ton-scale detector. 
With a 20 cm of Pb shielding backed up with 5 cm of pure cooper, the external gamma rate 
can be cut down to <10 Hz (>1 MeV threshold), a level that can be comfortably handled 
by commercial electronics. 



B. Time-correlated background rate 

As mentioned above, background in v-D charge current channel is hugely suppressed by 
the triple coincidence requirement. Here we focus on the time-correlated background in the 
IBD channel. We divide the background into two categories: reactor-associated and reactor- 
independent. They can be further sub-divided into correlated and accidental background. 



1. reactor-associated background 

On average, each fission produces 2-3 fission neutrons j^. Thus for a 60 MW th reactor, 
the total fission neutron flux from the core is estimate to be 5 x 10 18 /s. Most of reactor- 
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associated background is due to neutrons from the reactor core as well as secondary gammas 
from neutron capture on metal or concrete surrounding the core. Since CARR is a neutron 
scattering facility with many neutron guides from the core, it is difficult to accurately es- 
timate the neutron background at the detector location. As a start, we assumed an ideal 
spherical geometry and used a GEANT4-based toy Monte Carlo program to transport fission 
neutrons from the core through heavy water (1 m), water (1.65 m), and concrete shielding 
wall (2.1 m). The neutron spectrum emitting from the concrete wall is shown in Fig. [5] 
together with initial fission neutrons spectrum. The shielding factor for neutrons above 1 
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FIG. 5: Energy spectrum for the residual reactor fast neutron background predicted by a GEANT4 
Monte Carlo. 

MeV is calculated to be 3.8 x 10 19 under this geometry. About one third of these residual 
fast neutrons carry a kinetic energy >4 MeV (corresponding to a prompt visible energy > 
1 MeV). So without any neutron shielding around the detector, and taking into account the 
detector acceptance (~2x 10~ 3 at 7 m), we estimate a correlated rate of < 1 x 10 _4 /s, two 
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Pulse shape discrimination is reported 



order of magnitude lower than the IBD signal 
to be able to distinguish nuclear recoil from electron recoil signals in LS [30|, |3jJ . This could 
lead to another powerful background suppression factor. 

As mentioned earlier, fast neutron background is absent in the water and heavy water 
detectors. In all three detectors, slow neutrons will only make a single capture signal, con- 
tributing to accidental background only. Our high energy threshold on the delay-like events 
remove most gamma background from the natural radioactivity, but not single neutrons that 
captured in the target via n-Gd (LS or water) or n-Cl (heavy water) or high energy gamma 
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rays caused by neutron captured on envioronmental metal materials (Fe/Cr/Ni etc.) 32|. It 
is difficult to estimate delay-like back grou nd without direct reactor-on measurement. Just 
for reference, NUCIFIER experiment 



39[ estimated that such a background contribute to 



B/S of 1:1 for their IBDs [32 



with 10 cm of Pb shielding. A 20 cm Pb shielding at CARR 
may further suppress this background. 

It is interesting to note that other experiments at research reactors also have measured 
reactor-correlated background. For example, the ILL reactor neutrino experiment reported 
no such background [331 ] . 



2. reactor-independent background 

As the detector is placed at the surface, the reactor-independent background is domi- 
nated by neutrons created by cosmic muons (LS detector). Muons can be separated in two 
categories: a) "LS muons" with long trajectories in the LS so the detector itself has 100% 
tagging efficiecy, b) "corner muons" with short or no trajectory in the LS therefore missed 
by the detector. Fast neutrons produced by corner muons will be untagged, producing 
dangerous background to the experiment. 

To suppress untagged neutron background, we assume a simple outer muon veto system, 
e.g. two layers of plastic scintillator paddles, with an efficiency 95%. Let us further assume 
that the paddles cover an area of 4 m 2 on the surrounding dead materials or corners which 
would have been missed by a bare LS detector. The integral intensity of vertical muons at 
sea level is about 60m~ 2 s _1 sr^ 1 For an area of 4 m 2 , the muon rate is estimated to be 
60 x 4 = 240 Hz. If we veto all IBD-like candidate within 200 us to a muon detected by the 

n 

paddles [2J, additional deadtime introduced to the experiment is only 4.8%. 

The residual fast neutron background in the detector after the paddle veto is estimated 



using empirical parameterization in 34] 



N n = 4.14EJ' 74 x 10~ 6 /(/i ■ g/cm 2 ) . (8) 

where is the muon energy in GeV. Let us conservatively assume that the 5% unvetoed 
muons each have an effective path of 100 g/ cm 2 in the detector (although they only hit the 
dead surounding material), we get an untagged spallation neutron rate of 

240 x 5% x 4.14 x 4 a74 x 10" 6 x 100 ~ 0.01Hz 
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where surface muon average energy (E^) ~ 4 GeV []| has been been assumed. This back- 
ground is comparable to the IBD background. Other background due to muons hitting 
outside the 4 m 2 area is expected to be small. 

Cosmic ray induced fast neutrons in the LS can be further suppressed via following 
handles. First, neutron tagging efficiency can be improved by employing more layers of 
muon paddles. Second, as illustrated in [35| , a model independent way to remove all reactor- 
independent background is to use the "reactor on - off" . Third, as mentioned in Sec. IIVB 1\ 
pulse-shape discrimation technique may help to veto neutron recoil signals. For water and 
heavy water detectors, on the other hand, there is no correlated background of this nature. 

Muon spallation on carbon in the liquid scintillator will produce beta-delayed neutron 
emitters Li 9 or He 8 , two classical correlated background (/3=prompt, n=delayed) for the 
liquid scintillator-based experiments. Due to the long half-live (0.18 and 0.12 s), one can 
not tag them with the parent muons nor is it practical to veto all muons for a long ~s 
time window. If we extrapo.ate from the fitted results at Daya Bav near site: ~30/day/20- 
ton for an average muon energy of 55 GeV and rate of 20 Hz [36], we get a surface rate 
of 1.3/day/ton assuming that the production yield shares the same -E^' 74 dependence as 
the neutrons. This background is negligible compared to the IBD rate, and again, can be 
subtracted by reactor off data. 



V. SENSITIVITY TO STERILE NEUTRINO SEARCH 

In what follows, we investigate the sensitivity of these neutrino detector to potential 
anti-electron to sterile neutrino oscillations. For simplicity, we shall dwell in the 3 (ac- 
tive) +1 (sterile) framework. At short distance(<15 m) from the reactor, the oscillation from 
standard 3x3 mixing parameters can be neglected, so the electron anti- neutrino's survival 
probability is 

P ee = 1 - sin 2 sin 2 (1.27^Am 2 4 ) (9) 

in which L is the baseline (distance from reactor core to detector) in meter, E is the neutrino 
energy in MeV, and Am 2 4 is the mass square difference between V\ and z/ 4 mass eigenstates 
in eV 2 . Taking into account the oscillation, the detected neutrino spectrum for a perfect 



12 



detector (Eqn. HI [7J) is modified into 

N osc {E v ) = P ee N no _ osc {E u ) (10) 

In the reminder of this note, we chose a 60 MW^ reactor, 1 year running time, 1 ton 
fiducial mass, and 100% detection efficiency as the default exposure. Results for different 
exposure can be projected straightforwardly. 



A. Detector Response Function 

As mentioned in Sec. IIII Al and II II C[ the positron kinetic energy E e + is simply related 
to neutrino energy E v . On the other hand, the detectable prompt energy of the three media 
is different. For LS, as mentioned in Sec. IIII A[ we have E p = E e + + 1.022 MeV. For water 
and heavy water detector, due to the Cerenkov threshold, the annihilation energy is hardly 
visible, so E p = E e +. To get realistic visible energy spectrum, we convolve E p with a simple 
Gaussian smearing (resolution) and a step-wise threshold function, so 

N ms (E ms ) = T{E ms ) J N(E p )g(E ms - E p )dE p . (11) 



In this expression, Q(E V i S — E p ) is Gaussian with a width of 10%/ \fE p for liquid scintillator, 
and 30%/ \[E~ P for water and heavy water (Sec. IIII Bp . and 

T(E vis ) = 1 LS (12) 
= l{E vis > lMeV) 



0(E VIS < lMeV) 



water and heavy water . 



N(E P ) is the prompt energy spectrum assuming a perfect detector, i.e. iV„ _ osc in Eqn. H]or 
[7J for no oscillation hypothesis, or N osc in Eqn. [TU1 when disappearance is taken into account. 



B. x 2 Definition 

The sensitivity of an given experimental setup to a given set of oscillation parameter 
(sin 2 2# 14 , Am 2 4 ), in short, is the power that one could differentiate the measured spectrum 
from a non-oscillation spectrum. Typically one defines a \ 2 function as the measure of 
such difference. In this application, it should satisfy that 1) x 2 = when 6*14 = 0, 2) for 
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a given value of Am 2 4 , its variation w.r.t. sin 2 2^14 follows the standard \ 2 distribution. 
Then to determine the exclusion limit to a given confidence integral (e.g. 95.5%), for each 
given Amy we would scan over the value of sin 2 2# 14 to generate measured spectrum, and 
determine the boundary of the corresponding y 2 (e.g. \ 2 = 4). 



Omitting background related systematics [40j, the \ 2 can be defined as [37| 



* ^ AT* [1 4- rr 2 AT* I ^ ' 

i = l nbins vis,oscl ' b2b vis,osci 

+ (^-) 2 + (^-) 2 + ( — J' + (-) 2 • 

@norm @eshift &escale ®L 

where N % vis osc and N* is no _ osc , respectively, represent the ith energy bin in the visible energy 
spectrum (Eqn. [TTj) with and without oscillation. This is so-called "rate+shape" x 2 . The 
following systematics have been considered (see also Tabled]): 1) a 3% normalization uncer- 
tainty er norm (including reactor total neutrino flux, target protons, and detector efficiency) 
and its nuisance parameter a; 2) energy non-linearity including a shift cx es hift (0.02 MeV) 
and a scale factor a esc& \ e (1%), and their corresponding nuisance parameters r\ and f3; 3) 
2% bin-to-bin uncorrelated shape uncertainties <Jb2bi which is added to the denominator of 
the first term for simplicity instead of introducing Nu ns of pull terms; f l {f],P) represents 
fractional change of counts in bin % for a given set of parameter (rj,/3) away from (0,0); 4) 
a 10 cm position accuracy of the center of the core o^, conservatively being assumed to be 
along the radial direction, and the corresponding nuisance parameter 7. In principle the 
7 term can be absorbed in the normalization uncertainty, but we leave them here to ease 
later discussion on two-detector sensitivity. The effects of detector resolution and threshold 
have been included automatically by using detected N* is osc (E vis ) and N l vis non _ osc (E vis ) (see 
Eqn. [T2l) in Eqn. [13j If one wants to perform a "rate-only" analysis, it is equivalent to using 
the above x 2 with a single visible energy bin and set rj, (3, and a^b to zero. 

The impact of reactor flux can be further suppressed if we choose to use two identical 
detectors located at two different baselines, similar to the set up in the Daya Bay and RENO 
experiments 

aa 

. In this case x 2 can be redefined as 

[Niiosc - N^no-osc -(l + a + e d + + f(r) d , (3 d ))f 

d=n,f;i=l,nbin 



2 ^ — v L vis,ase vis,no—osc V 1 1 1 (L a )' 2 J \« //J , A .» 

X ~ ^ [1 + rr 2 N d,i 1 

vis,osc[ 1 w b2b vis,osci 



« ^2 , , 7 ^ 2 , v- r/ e d , 2 , , rj d N2 , , /3° 



+(-^) 2 +(-^) 2 +$:[(-^) 2 + (-^) 2 + 

'-'norm ®h . <-f f Oeshift Oescale 
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^ 2 1 ; 



in which the superscript d runs between "near" and "far" to represent different quantities 
for the two detectors. We have also added a detector uncorrelated efficiency uncertainty a d ^ 
(0.5%) and its corresponding nuisance parameter e d . A summary of systematic components, 
the values, as well as whether they are correlated between the two detectors is given in 
Table E 



Systematic uncertainty 


Value assumed 


Nuisance parameter 


Near-Far Correlated 


Overall normalization a norm 


3% 


a 


C 


Detector relative efficiency & e ff 


0.5% 


e d 


u 


Energy shift a es hift 


0.02 MeV 


rf 


u 


Energy scale a esca i e 


1% 


f3 d 


u 


Reactor spectrum shape ovib 


2% 




u 


Baseline ai 


10 cm 


7 


c 



TABLE I: Summary of systematic effects included in the two-detector x 2 function in Eqn. [TH The 
last column indicates whether the component is correlated between the detectors: C=correlated, 
U=uncorrelated. Same systematic uncertainties are assumed for single-detector x 2 (Eqn.fTB"]) except 
for detector efficiency, which we have absorbed into the normalization uncertainty. 

One should note that a couple of conservative approximations have been made in Eqn. [TH 
First, instead of introducing nbins nuisance parameter for the bin-to-bin shape uncertainties, 
we assumed that these uncertainties are also uncorrelated between the near and far detectors, 
and lump them to the denominator of the first term just like the statistical uncertainties. 
Second, we have omitted detector correlated energy shift and stretch, and have assume 
rj d , (3 d as detector uncorrelated nonlinearity. Both approximations have been verified to 
have negligible impact to the sensitivity results. 

C. Baseline optimization 

Reactor neutrinos have a energy spectrum ranged up to 9 MeV, as shown in Fig. [TJ The 
IBD (v-D) neutrinos has a energy peak at E v ~ 3.7(5.7) MeV. If the true Amy is around 
1 eV 2 , as implied by the global analysis, naively one would put the detector close to the first 
oscillation maximum to maximize the analyzing power, i.e. L osc = f 1 yfj^m?. > translating 
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into 4.6 (7.0) m for IBD and v-D CC neutrinos. 

However, the above discussion is incomplete since we have omitted influence from statis- 
tics. The fact that the event rate is inversely proportional to L 2 makes the optimal baseline 
deviates from the naive oscillation maximum. 

A more elaborated analysis was made by employing \ 2 definition from Sec. IV Bl We 
assumed a fix parameter pair (sin 2 (2# 14 ) = 0.1 and Am 2 4 = 1 eV 2 ), and our later conclusion 
does not change significantly with the value of sin 2 (2^14). Energy thresholds in Eqn. [12] have 
been assumed, but for simplicity we assumed no energy smearing and set nuisance parameters 
for energy nonlinearity (77, /3), baseline uncertainty (7), and the bin-to-bin uncertainty a^b 
in Eqn. H31 to zeros. The optimal baseline was determined by scanning through baseline 
to find the maximum of x 2 - 111 Fig- El the value of \ 2 vs - baseline is shown for all three 
type of detectors with a "rate-only" or "rate+shape" analysis. Coincidentally, we observe 
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FIG. 6: Baseline optimization for single detector in "rate only" and "rate+shape" analyses for LS, 
light water and heavy water, with Am 2 4 = 1 eV 2 and sin 2 2#i4 = 0.1. See text for details. 

that 7 meter is sufficiently close to the best baseline for all three detection methods either 
in "rate+shape" or "rate only" analysis. \ 2 map for different values of baseline and Am 2 4 
have also been shown in Fig. [71 One observes a general trend that the sensitivity for larger 
value of Am 2 4 increases with decreasing far detector baseline, which is expected from the 
LAm 2 /E dependence of the shape. 

For reference, number of detected neutrino events at CARR using three different targets 
(1 ton-1 year-100% efficiency) at 7 m are tabulated in Table HTl 
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Sensitivity for scintillator with single detector ~| Sensitivity tor light water with single detector 




Sensitivity for heavy water with single detector | 




baseline(cm) 



FIG. 7: x 2 ma P vs - baseline and Araj 4 for LS (top left), light water (top right) and heavy water 
(bottom) in "rate+shape" analysis. 3% normalization uncertainties has been considered in the x 2 
calculation, but not other systematic uncertainties. 

Target Liquid scintillator H2O D2O 

Non-osc (/ton/year) 420129 369505 5515 

Osc (/ton/year) 394667 345329 4998 

TABLE II: Event rate at 7 m with 100% efficiency. Non-osc: sin 2 20 14 =0; Osc: Am 2 4 = 1 eV 2 , 
sin 2 26*14=0.1. 

For the two-detector scenario, since 7 m is approximately the closest distance that we 
can put the detector, we settle the near detector at this baseline. To determine the optimal 
baseline for the far detector, instead of using the two-detector \ 2 in Eqn. [TH (in which 
one still have to input the expected "non-osc" reactor spectrum), we adopted an approach 
to construct a pure relative measurement. We used the one-detector x 2 111 Eqn. [131 an d 
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demanded that N l vi osc 



(L n ) 2 n>i 

\2 vis,osci 



a scaled near detector spectrum. In general, this 



approach does not give the best sensitivity as it entirely omits the theoretical knowledge on 
reactor neutrino spectrum. On the other hand, the systematic uncertainty due to theoretical 
assumption is also completely avoided. The results of the baseline scan under this approach 
are shown in Fig. [8]). The best far detector baseline is about 11 m for both LS and water 
in "rate+shape" analysis, and about 9 m for LS and 10 m for water in "rate only" analysis. 
The optimal baseline for heavy water is about 14 m for both analyses. For two-detector 
discussions in the rest of this paper, we will assume the optimal baselines in "rate+shape" 
analyses, i.e. a near detector at 7 m, and a far detector at 11 m (LS & water) or 14 m 
(heavy water) [41 ]. 

baseline optimization for the far detector 
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FIG. 8: Baseline optimization for far detector with the near one placed at 7 m from the core in 
"rate only" and "rate+shape" analyses for LS, light water and heavy water, with Am 2 4 = 1 eV 2 
and sin 2 2# 14 = 0.1. See text for details. 



D. Sensitivity results 

1. Single detector 

Using the full x 2 definition in Eqn. [131 t ne 95.5% exclusion curves for the three detection 
methods are shown in Fig. [9] in the (sin 2 2# 14 , Am 2 4 ) plane. 

With "rate-only" analysis, we observe that for all three techniques, the 2a sensitivity 
limits are around 0.1 even at optimum Am 2 4 , due to the large normalization uncertainty. 
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FIG. 9: Sensitivities with a single detector using LS, water and heavy water in "rate+shape" 
analysis. All systematics in Eqn. [13] have been considered. 

The sensitivity curves for LS and light water track with each other, although the latter is 
slightly worse due to a larger detection threshold. Both experiments are most sensitive to 
Am 2 4 ~ 0.65 eV 2 . For Am 2 4 larger than ~ 3 eV 2 , the oscillation as a function of energy 
becomes so fast that it get smeared out by the energy resolution of the detector. In this 
case, one measures a constant deficit (oc 1/2 sin 2 (2# 14 )) independent of the baseline, giving 
rise to a constant sensitivity at large Am 2 4 . The heavy water sensitivity is not so much 
worse than the other two, as the dominating uncertainty comes from the normalization, not 
the statistics. The value of Am 2 4 where the best sensitivity occurs is ~1 eV 2 , higher than 
that of the IBD, due to higher reaction threshold (~4 MeV) thereby higher average detected 
neutrino energy. 

The situation is drastically improved with a "rate+shape" analysis. All three exclusion 
curves moved to much smaller value of sin 2 2# 14 in Fig. [9J LS and water in particular, as the 
constraints from the shape will seriously combat the large normalization uncertainty. It is 
also interesting to note that the value of Am 2 4 where experiments are most sensitive to has 
undergone significant changes compared to that in "rate-only" analysis. This can also be 
understood as an effect from the shape constraints. For example for IBD, Am 2 4 ~ 0.65 eV 2 
the overall disappearance in rate is the largest, on top of the 3% uncertainty in normalization. 
However, at Am 2 4 = 0.65 eV 2 the shape distortion is flatter compared to say Am 2 4 = 1, 
so it would be relatively easier to choose a normalization nuissance parameter to balance 
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spectrum distortion. For heavy water, the "rate+shape" analysis helps, but not as much as 
the IBD, as a result of lower statistics in each energy bin. 

2. Two detectors 

The sensitivity with two detectors can be investigated in a similar way using the \ 2 m 
Eqn. HU Conceptually, unlike the design of Daya Bay and RENO experiments p, [J (using 
well-known Am 2 ), with an unknown Am 2 4 the "rate-only" relative measurement becomes 
much dicier. For certain values of Am 2 , the normalized event rates at near and far sites 
would equal to each other therefore would cancel the sensitivity in the near/far ratio. Under 
this "unlucky" situation, one could still gain some sensitivity back by relying on the flux 
prediction, but the main purpose of two-detector design would be largely undermined. 

The story is drastically different once the detector shape information is used, due to the 
L/E dependence of oscillated spectrum. Within two-detector scenario, the sensitivity curves 
with "rate+shape" analysis are shown in Fig. (TUlfor LS, H 2 and D 2 0. One clearly observes 
an improved sensitivity compared to Fig. [9] for all detectors. This is anticipated, since the 
X 2 construction in Eqn. [TH has used information from both detectors as well as the flux 
prediction. 

Sensitivities with two detectors 

> 

% 

to 

10 



10" 1 

sin 2 29 

FIG. 10: Sensitivities with two identical detectors using LS, water and heavy water in "rate+shape" 
analysis. All systematics in Eqn. [T5] have been considered. 

To study contribution from each individual systematic component, we compared the 
variation of the exclusion curve on (sin 2 29^, Am 14 ) when "turning off" systematics one- 




20 



by-one. The results for the LS detector (most sensitive one) at Am 2 = 1 eV 2 is shown in 
Table Hill Separated systematic uncertainties by taking the quadrature difference [42 1 are 
also tabulated in the table. 



Errors considered 


Single detector 


Two-detector 


A+B+C+D+E+Stat. 


0.02748 


0.02017 


A+B+C+D+Stat. 


0.02701 


0.01967 


A+B+C+Stat. 


0.02680 


0.01934 


A+B+Stat. 




0.01746 


A+Stat. 


0.02321 


0.01691 


Stat, only 


0.01055 


0.00885 



Stat Normalization Efficiency E stretch E shift Baseline 
Single detector 0.011 0.021 - 0.013 0.004 0.005 

Two-detector 0.009 0.014 0.004 0.008 0.004 0.004 

TABLE III: Upper table: 2-cr sensitivity on sin 2 20±4 when adding in systematic effects one-by-one. 
Systematic uncertainties include: (A) normalization, (B) detection efficiency, (C) energy scale 
stretch, (D) energy scale shift, and (E) baseline uncertainty. The statistical and bin-to-bin shape 
uncertainties are always combined into a "stat" uncertainty. Target=LS. Am 2 = 1 eV 2 . Single 
detector: 7 m. Two-detector: near i 7 m, far @ 11 m. Lower table: breakdown of uncertainties 
by taking the quadrature differences from the upper table. The efficiency uncertainty for a single 
detector has been included into the normalization uncertainty. 

One sees that two-detector scheme not only improves the statistics, but also help to re- 
duce the normalization uncertainty. The fact that the normalization uncertainty does not 
disappear completely is a result of the interplay between the rate and shape constraints. 
From experimental point of view, if detectors are constructed as movable, one could con- 
sider a "swap" between the near and far detector in order to further suppress systematic 
uncertainties. 
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VI. DISCUSSION AND CONCLUSION 



Under the framework of "3-active+l-sterile" neutrino mixing, we have conducted a study 
of the sensitivities to sin 2 2# 14 at short baseline (< 15m) to a research reactor (CARR) using 
three targets (LS, H 2 0, and D 2 0). This study suggests that in the absence of background, 
the LS detector has the best sensitivity due to higher IBD reaction rate and superior energy 
resolution. For experiment detecting v— D CC scattering using heavy water, the event rate is 
far less, with a huge suppression of potential background nevertheless. From the comparison 
between "rate-only" and "rate+shape" analyses, we conclude that the spectrum distortion 
provide crucial handle to oscillation therefore a key to the experiment design. To cancel 
uncertainties from the reactor flux prediction, we compared the performance of a single- 
detector and two-detector design. The latter leads to not only a better sensitivity but also 
a suppression of systematic uncertainty. Under the current best scenario (liquid scintillator, 
two-detector, no background, and "rate+shape" analysis), a ton-scale detector can reach a 
sensitivity (95.5%) of sin 2 26* 14 of order 0.02 over a range of Am 2 4 between a few tenth to a 
few eV 2 , which covers a 5a discovery range to the value suggested by global analysis. 
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